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Abstract. The interaction between a high-frequency dilational mode of a thin dielectric film 
and an optical cavity field is studied theoretically in the membrane-in-the-middle setup. A 
derivation from first principles leads to a multi-mode optomechanical Hamiltonian where 
multiple cavity modes are coupled by the thickness variation of the membrane. For membrane 
thicknesses on the order of one micron, the frequency of this dilational mode is in the 
GHz range. This can be matched to the free spectral range of the optical cavity, such 
that the mechanical oscillator will resonantly couple cavity modes at different frequencies. 
Furthermore, such a large mechanical frequency also means that the quantum ground state of 
motion can be reached with conventional refrigeration techniques. Estimation of the coupling 
strength with realistic parameters suggests that optomechanical effects can be observable 
with this dilational mode. It is shown how this system can be used as a quantum limited 
optical amplifier. The dilational motion can also lead to quantum correlations between cavity 
modes at different frequencies, which is quantified with an experimentally accessible two- 
mode squeezing spectrum. Finally, an explicit signature of radiation pressure shot noise in this 
system is identified. 



PACS numbers: 42.50.-p, 42.50.Wk, 78.20.hb, 42.65.Yj 
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1. Introduction 

Research on mechanical systems in the quantum regime has reached several milestones in the 
past few years. A micromechanical oscillator in the quantum ground state was first observed 
in an experiment by O'Connell et al. HI, where the high-frequency dilational motion of a slab 
of dielectric material was coupled to a superconducting qubit. This enabled verifying that 
the mechanical oscillator was in the ground state, as well as the ability to controllably create 
single phonons. The very high mechanical resonance frequency of 6 GHz ensured that the 
quantum ground state was reached simply by conventional refrigeration techniques. 

Quantum behaviour of mechanical oscillators has later been achieved in cavity opto- or 
electromechanical systems, where the mechanical motion is coupled to light or microwaves. 
Resolved sideband cooling of mechanical motion to the quantum ground state [[2], [31 was 
achieved in an experiment by Teufel et al. on an electromechanical system Hi. This technique 
has also been successfully employed in the optical regime with silicon nanoscale resonators 
[El, where the quantization of mechanical motion has been directly measured (6). In a 
setup coupling light to the collective motion of cold atoms, quantization of motion [7] and 
ponderomotive squeezing of the light [8J have been reported. 

The optomechanical system of a thin silicon nitride membrane placed inside an optical 
Fabry-Perot cavity [0 has also attracted considerable attention recently IfTOl ITTl [T2l [T3l [T4l 
[T51 QSD. The flexural modes of the membrane give rise to the same physics as a cavity with 
a movable end mirror. In addition, one can achieve coupling of the light to the square of 
the mechanical displacement flUO, which could enable nondemolition measurements of the 
phonon number [TTOll or of phonon shot noise [fT2l . The idea of using such a membrane as a 
transducer between optical and electrical degrees of freedom has been proposed [16]. Finally, 
the membrane's coupling to light was recently exploited in an experiment to effectively couple 
the motion of the membrane and the motion of cold atoms IfPTll . 

While the membrane-in-the-middle setup is promising, laser cooling of the flexural 
modes to the ground state is a considerable challenge. One reason for this is the typically low 
resonance frequencies of these modes, usually in the range of 0.1-10 MHz, which corresponds 
to a very low effective temperature necessary to reach the ground state. However, in addition 
to the flexural modes, the membrane also experiences thickness oscillations, i.e. a dilational 
mode similar to the motion studied in the experiment by O'Connell et al. [1]. If the membrane 
is very thin, the thickness oscillation has a much higher frequency than the flexural modes. 
The frequency is approximately given by v B /(2do), where v s is the speed of sound in the 
membrane and d is the membrane thickness. Note that thin-film dilational modes of this type 
have been detected in free-standing silicon membranes by using laser pulses |[T8l [T9ll . For 
silicon nitride membranes with thicknesses of a micrometer or less, the mechanical resonance 
frequencies will be on the order of GHz or tens of GHz. Hence, provided that the coupling 
between this dilational mode and light is sufficiently strong, the membrane-in-the-middle 
setup could also give rise to an optomechanical system where the mechanical oscillator is 
in the ground state already at dilution refrigeration temperatures. Another useful feature of 
the very high mechanical frequency is that it is of the same size as the free spectral range 
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in standard optical cavities, so the dilational mode will effectively couple cavity modes at 
different frequencies. This type of frequency matching was recently exploited in a different 
type of optomechanics experiment [|20l . 

In this article we investigate the coupling between the thickness variation of a thin 
dielectric membrane and an optical cavity field. In Sections [2] and |3} we characterize the 
uncoupled optical and mechanical modes, respectively. The coupling comes about both due 
to the membrane thickness variation itself (a boundary effect) as well as the fact that elastic 
strain in the membrane leads to a change in the electric permittivity (a bulk effect). In general, 
both effects must be taken into account when light interacts with the motion of dielectric 
objects f20l I2TI . In Section |4j we derive a second quantized Hamiltonian and quantify the 
coupling strengths between the optical cavity modes and the mechanical oscillation. By 
inserting realistic parameters, we show in Section [5] that the coupling can be strong enough 
for optomechanical effects to be observable, provided that the as yet unknown mechanical 
quality factor exceeds a certain value. In Section [6j we show how this system can be used 
as a quantum limited optical amplifier, due to the vanishingly small thermal occupation 
of the mechanical mode. Finally we study how one can produce quantum correlations 
between different cavity modes with this setup in Section [7] We calculate a two-mode 
squeezing spectrum and find that one can achieve considerable squeezing for sufficiently large 
mechanical quality factors. In addition, we identify an explicit signature of radiation pressure 
shot noise. 

2. Optical cavity modes 

We will consider the setup of a thin dielectric membrane placed inside an optical Fabry-Perot 
cavity, depicted in Fig. [T] The length of the cavity will be denoted by L, and the equilibrium 
membrane thickness by do. The center of the cavity is placed at x = 0, and the right (left) 
mirror is positioned at x = ±L/2. The electric field inside the cavity can be expressed in 
terms of orthonormal modes <fr/(r), 




where ui is the mode frequency, £ is the electric permittivity in vacuum, and a\ are photon 
annihilation (creation) operators with [a;,a|,] = 5ij'. In the absence of the membrane, the 
functions 3>j(r) are the Hermite-Gauss modes [1221 . In general, the subscript / represents a set 
of quantum numbers. However, we will restrict ourselves to the Gaussian TEM o-modes only 
and consider polarized light. In that case, I = 1,2, ... will be the longitudinal mode index. 

The width of the cavity modes changes on the length scale of the Rayleigh range 
xr = kiwl/2, where ki is the wavenumber and w the waist size of the cavity mode. We 
will concentrate on the situation where the membrane is positioned at xm well within the 
Rayleigh range from the waist, such that |xm| <C £r. In this case, the mode functions inside 
the membrane can be approximated by <&z(r) = $/(r)v [23 ] where v is a polarization vector 
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Figure 1. Proposed experimental setup, a) A thin dielectric membrane inside an optical Fabry- 
Perot cavity. The membrane is placed close to the waist of the cavity mode, b) Close-up view 
of the membrane when the membrane thickness is at its equilibrium value do. The waist radius 
of the cavity mode is wo and the side length of the square membrane is b. c) Vibrational modes 
in the membrane leads to thickness fluctuations around the equilibrium value do- 
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We have introduced ry = (y, z) as the component of the position vector in the yz-plane, and 
Um is the membrane's index of refraction. The wave number is ki = un/c, with c being the 
speed of light. The dimensionless coefficients C\ and the phases (pi follow from normalization 
and boundary conditions. We do not include their expressions here, but note that C\ is of order 
1 and 0;/ — 0z ~ {V — l)n/2 if the membrane's intensity reflection coefficient 
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is much smaller than 1. For Si 3 N 4 , n M 2, such that R{k) oscillates between and 0.36 as a 
function of kdo. 

The mode frequencies are given by coi = cof + 5u)i, where uf = lirc/L are the evenly 
spaced frequencies in absence of the membrane for large longitudinal mode numbers /. The 
correction to the frequencies in presence of the membrane can be expressed in terms of 
k® = u}f/c and becomes E31 
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The frequencies have an oscillatory dependence on the position of the membrane [|9l whose 
magnitude depends on the reflection coefficient R(k). Note also that the shifts for neighboring 
longitudinal numbers I and I + 1 have opposite sign. The contribution /3(kf ) has only a weak 
dependence on I and can be thought of as an unimportant shift of all frequencies. The mode 
frequencies as a function of membrane position are shown in Fig. |2j 
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Figure 2. Cavity mode resonance frequencies in units of wfsr = cir/L as a function of 
membrane position %. We have used a membrane thickness of do = 1 /im, a wavelength of 
A = 1064 nm, and an index of refraction = 2, giving R(k) = 0.2. The frequencies have 
an oscillatory dependence on xm, which means that they are in general not evenly spaced. The 
blue dashed line to the right is placed at a membrane position where the frequency spacings 
Si = u)i — differ, i.e. Si ^ Si + i ^ <5;+2- The red dashed line to the left is placed at a 
high-symmetry point where Si = Si + \ = 8i +2 - 



3. Dilational modes of the membrane 

In this section, we study the high-frequency dilational modes of the membrane, which is 
assumed to be square with side lengths b and an equilibrium thickness d . The membrane lies 
in the yz-plane with its front and back surfaces at x± = xm ± d /2. The displacement field 
in the membrane will be denoted by u(r). Assuming that the membrane material is isotropic, 
the displacement eigenmodes are given by the equation 



/jV 2 u + (K + | J V ( V ■ u) = -p lo 2 u, (5) 

where K is the bulk modulus, \i is the shear modulus, p is the mass density, and u is the mode 
frequency. To find the displacement field, we also need to specify the boundary conditions, 
which depend on the details of how the membrane is connected to its mechanical support. For 
simplicity, we will work with the so-called simply supported boundary conditions 

u x = u y = cr zz = , z = ±6/2 (6) 
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where an are the diagonal elements of the stress tensor. This means that the transverse 
displacement as well as the normal component of the stress on all surfaces are zero. This 
simplification might not be the most realistic choice of boundary conditions. However, we do 
not expect that other boundary conditions will bring significant changes to the properties we 
are interested in, namely the frequency spacings between the modes and the optomechanical 
couplings. 

The boundary conditions and the equation of motion are satisfied by the field u(r) oc 
iimnfr) with 



u, 



U, 



K x sin 



K y cos 



K~ cos 



[x 



"0 
71 

d 

7T 



y X - X M 



x 



Xm) 



sm 



cos 



sm 








sin 





(7) 







(■* 




sin 




















cos 













where Kj are constants, m, n 



1,2,3..., and we have restricted ourselves to the lowest 

is given by 



non-zero wave number (ti /do) in the x-direction. The angular frequency u = oo mn 



(8) 



(K + 4/i/3)/p for longitudinal waves and 



where the propagation velocity v — v s = 

v = y/jjp for transversal waves. 

Since the thickness do is much smaller than the width b, the longitudinal vibrational 
modes will consist mostly of displacements in the x-direction, whereas the transverse modes 
consist mostly of displacements in the plane of the membrane. In the remainder of this 
article, we will focus on longitudinal modes only. We emphasize that the transverse modes 
can also couple to the optical cavity modes through the photoelastic effect. However, since 
the longitudinal and transverse mode frequencies do not coincide, we can safely ignore the 
latter below. The normalized longitudinal modes are given by Eq. ^ with K y /K x = mdo/b, 
K g / K„. = ndo/b, and K x = (1/(1 + (m 2 + n 2 )(d /b) 2 )) 1/2 w 1. Note that for these modes, 
both the normal and shear stress on the surfaces x = xm ± d /2 are zero, as they should be 
on a free surface. 

Ignoring the transverse modes, the total displacement field becomes u(r) = 
J2mn l )mnUmn( r )' In other words, u(r) is the displacement field of the vibrational modes 
with angular frequencies close to irv B /do. The length r] mn characterizes the amplitude of the 
mn-mode. We can quantize the displacement field by letting r] mn — >■ r/ zp f irrm (c mn + c^J, 
where c mn is a phonon annihilation operator for the mn-mode. The zero point fluctuations 
are characterized by the length 77 zp f i?Tm = \J h / (2muj mn ) , where m = pd b 2 /8 is the effective 
mass of each mode. 
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We would like to study the coupling between the optical cavity modes and a single 
vibrational mode, namely the fundamental m — n — 1 mode. This is possible if the modes are 
well separated in frequency, which requires that the mechanical quality factor Q ^> (b/do) 2 . 
This is one of several reasons for using a membrane with a small area to reduce the aspect 
ratio b/d . If this requirement for the mechanical Q is too strict, it is however also possible 
to focus solely on the fundamental mode if the light couples much more strongly to it than to 
the other longitudinal modes. We will address this issue in the next section. 



4. The Hamiltonian for the coupled system 

We now derive the Hamiltonian for the coupled optomechanical system. We assume that we 
can focus on a single mechanical mode with m — n — 1. We let uu — > com ~ v s 7i/do, 
Vmn — > V' an d ~~ * u, such that the displacement field is u(r) = r^u(r). We define 
ri = (rr_, rii) and r r = (x + , rn), which describe the equilibrium left and right membrane 
surfaces, respectively. Note that the membrane thickness at in-plane position ry is given by 
d(rn) = d + 2u x (r r ). The displacement field will later be quantized as described in Section 
[3] The coupling to other mechanical modes will be commented on in Section 4.3 



4.1. Derivation of the full Hamiltonian 

For a static membrane, the Hamiltonian is H = J2i fujJia\ai. In the presence of membrane 
vibrations, the Hamiltonian becomes H = H + Hi + H 2 + Hm, where is the energy 
of the mechanical fluctuations. The term Hi is the energy change of the electromagnetic 
field due to the change in membrane thickness, whereas H 2 is the energy change due to the 
photoelastic effect, i.e. the fact that a strain in the membrane gives rise to a change in the 
electric permittivity tensor. Using perturbation theory ll25l . these terms can be expressed as 
Hi = 1/2 / drii/ij(r||) where 

&i(r||) =t7u a ,(r r )(e M , -eo) [E 2 (n, t) + E 2 (r r , t)l (9) 



E (r,t) . 

r?=0 



Note that the expression for the surface contribution Hi is only well-defined when the electric 
field is continuous across the boundary [25] and the optical modes are non-degenerate, both 
of which are satisfied here. The strain-dependent permittivity in the membrane is denoted by 
£m(i~), and Em,o = ^oy^M ^ s i ts equilibrium value. We have assumed that the amplitude t] 
is small compared to the wavelength of the light and the equilibrium thickness do- We also 
assume that the optical cavity modes are polarized in the yz-plane, such that £m(t) is the 
in-plane membrane permittivity. Note that when c^mM / 'dr)\ v=0 < 0, i.e. when the electric 
permittivity decreases as the thickness increases and vice versa, the contributions Hi and H 2 
have opposite signs. 

In the presence of strain described by the strain tensor e^, the deviation of the 
inverse electric permittivity tensor e^ 1 from its equilibrium value is given by 5(e^ 1 )ij = 
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Y,kiPijkieki/£o, where pijki is the photoelastic tensor [|24]|. Assuming that the membrane 
is isotropic in the yz-plane and that x,y,z point along the principal axes of the membrane 
material, we have (e^ 1 )^ = s^/Sij with e yy = e zz = Sm- Neglecting the small displacement 
components u y and u z , we then arrive at deui^) / d7]\ v =o ~ — e^p/so du x (r)/dx, where 
p = Pyy XX and u x (r) is given by Eq. |7]). For Si0 2 , p yyX x = 0.271 11241 . We are not aware 
of any measurements of the photoelastic coefficients for Si3N 4 , but it seems reasonable to 
assume that they are not very different for Si 3 N 4 than for Si0 2 . In any case, they are expected 
to be of the same order of magnitude. 

We now quantize the displacement field by letting r] — > r] zpi (c + cM , where c and & are 



phonon annihilation and creation operators and ?7 zp f = ^h/(2rhuou) is the size of the zero 
point fluctuations. By inserting Q and |7]) and ignoring terms that do not conserve photon 
number, we arrive at the multi-mode optomechanical Hamiltonian 

H = ^2 ftuia\cn + fiouuc'c + KTiji {c + c f ) a\av . (10) 
i i,v 

Below, we will focus on the intermode coupling = r i+lji between adjacent longitudinal 

optical modes. Defining g\ = coir] zp f/L and inserting the expressions for Ci and <pu we find 

T M+1 = giAu +1 , where 




. , 1 \ ^(fcQcot(n M Mo) , 

AM+1 = - 1 1 " l-(2n M /,y) (l-^ ) cos 2 (2^ M)) ^ Sm(2 ^ M) ' (U) 



The quantity is a transverse overlap integral which is defined in Eq. (13) below. Note 



that sj R{ki) cot(n M kid ) — > when cot(n M ^o) — > ±oo, and that cot(n M /c;<io)/(l — 
(2nM^c?o/ 7r ) 2 ) - > 7r /4 when nyihdo — >■ tt/2. Thus, except for special values of the 
membrane position x M or the thickness d for which the coupling vanishes, the dimensionless 
number A i ;+1 is always of order 1. The coupling strength is therefore determined by gi, 
which is the familiar coupling rate in the canonical optomechanical system of a cavity with 
a movable mirror. This also determines the size of the intramode coupling 1^. We have 
used ~ h+i^M, yj^i^i+i ~ and — 0/ tc/2 to simplify the expression for 
the coupling strength. These are very good approximations for large longitudinal numbers I 
and small reflectivities R(k{). For further details on the derivation of Eq. (11), we refer to 
Ref. 11231 . where the calculation of coupling strengths to the flexural modes produce similar 
expressions. 

The first term in the paranthesis in Eq. (11) comes from the thickness variation whereas 
the second from the photoelastic effect. We see that these terms are of opposite signs for 
do < ir/(2nM.ki) an d °f me same sign for do > Tc/(2riMki). This does not contradict the 
fact that Hi and H 2 always have opposite signs for p > 0, since we are now considering 
the coupling between two particular optical modes, and not the entire Hamiltonian. We also 
see that the two terms are of comparable size for reasonable parameters. Finally, note that 
the photoelastic term vanishes when the membrane thickness becomes much larger than the 
optical wavelength. This is as expected, since this term is, roughly speaking, an integral of 
two orthogonal optical modes over the thickness of the membrane. 
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4.2. Reduction to two optical modes 

When focusing on a single optical mode with index / and the intramode coupling only, 
the Hamiltonian is reduced to that of a standard optomechanical system. However, for thin 
membranes, the frequency of the dilational mode u>m will typically be very large compared 
to a typical cavity linewidth, denoted below by k. Indeed, the mechanical frequency can be 
matched to the free spectral range of the cavity. Hence, we will consider two optical modes / 
and / + 1 that are separated by the mechanical frequency wm, such that the intermode terms 
proportional to r^ +1 and r i+lj / are resonant. This means that the annihilation (creation) 
of a phonon will convert a photon from mode I (I + 1) to mode / + 1 (I). Defining 
5 1 = uoi — note that the frequency shift ([4]) ensures that <^ +1 ^ 8\ « <5/ +2 whenever 
cos(2fc;a;M) ~ cos(2/c°xm) ^ (see Fig. [2]). Thus, by placing the membrane at a position 
where \6i+i — 8\\ > k, we can focus on just two optical modes. Furthermore, in the resolved 
sideband regime um k, we can neglect the intramode coupling as well as the counter- 
rotating terms in the intermode coupling. 

The maximal difference between <^ +1 and 8i occurs when cos(2/^xm) = ±1- However, 



according to Eq. (11), the coupling rate r^ + i is zero at these points. This is not expected to 
cause a problem in practice since maximizing — 5i\ is not really necessary. One only 
needs to ensure that |<^ +1 — 8i\ 3> k. To show that this is feasible, let us for example assume 



that cos(2kfxu) = 1/2 and \JR{k) = 0.2. By expressing the linewidth k = en /(LF) by the 
cavity finesse J 7 , we get \8i +1 — 8i\/k ~ 0.1 J 7 3> 1 for any reasonable finesse. Hence, we 
can safely focus on two cavity modes and neglect the others. The vertical blue (dashed) line 
in Fig. [2] is placed at a membrane position where both 5i + i ^ 5i, 5i +2 and the coupling rate 
T^+i is nonzero. 

For a given I, we let a/ — > Si, — > a.2, and = T, which leads to a model 

where two longitudinal cavity modes are coupled by the fundamental dilational mode of the 
membrane, 

2 

H = ^2 ^Wjajdj + huyi&c + hT {c d^ai + c^o{a 2 ) , (12) 

i=l 

with uj 2 — u)\ ~ t^M- The indices 1 and 2 are just for labeling the cavity modes and do 



not represent the actual longitudinal numbers 1,1 + 1 ^ 1. We note that the model (12) 
was recently studied in the strong coupling regime r ~ k ll26l l2~71l . where e.g. quantum 
nondemolition measurements of photon and phonon numbers become possible. The coupling 
in the experimental setup we study here is however not strong enough to see these effects. 

4.3. Coupling to other mechanical modes 

In the derivation of the interaction Hamiltonian we ignored all mechanical modes except the 
longitudinal fundamental mode with m — n — 1. However, as noted in Section [3j unless 
the mechanical quality factor Q ^> (b/do) 2 , this can only be justified if the optomechanical 



coupling to other mechanical modes close in frequency is significantly smaller. In Eq. (11), 



we expressed the coupling rate to the fundamental mode by the overlap integral Ji i, defined 
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where the integral is over the transverse area of the membrane and we have assumed that the 
cavity waist is centered on the membrane. The integral In is always of order 1. For other 
longitudinal mechanical modes that are close to the fundamental in frequency, i.e. modes 
with arbitrary but moderate values of m and n, the coupling will be given by the same 
expression as in Eq. (fTT|), except for the substitution In — > I m , n , with 

71 ( b\ 71 



71 Wq J A 
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b V ' 2 , 

Symmetry dictates that the coupling to modes with m or n even becomes zero. In Fig. [3j we 
plot the overlap integral I m>n normalized to the value of I\ t \ for the three modes closest in 
frequency to the fundamental mode as a function of cavity waist size. We see that when the 




2w /b 



Figure 3. The transverse overlap integral I mn normalized to the integral for the fundamental 
mode /i i as a function of the cavity waist diameter 2w$. This shows that when the waist size 
becomes comparable to the membrane's transverse dimension b, coupling of the light to the 
fundamental mode m = n = 1 is significantly stronger than to the other vibrational modes. 



cavity waist is comparable in size to the membrane's transverse dimensions, i.e. 2wq ~ b, the 
optomechanical coupling to the fundamental mode is dominant, which is another reason for 
using a membrane with a small side length b. The neglect of the other modes with morn > 1 
thus appears to be justifiable, although one should keep in mind that they are there and might 
influence experiments. 



4.4. Strong optical drive 

In Section [5j we will show that the coupling strength Y in Eq. (12) is typically very small 
compared to the cavity and oscillator decay rates. We therefore consider the situation where 
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one of the cavity modes is strongly driven. Let us first assume that the upper mode (mode 2) is 
strongly driven by a laser with frequency w D ~ cu 2 . Transforming to the frame rotating at the 
drive frequency wd, and neglecting small nonlinear terms, we get the parametric amplification 
Hamiltonian 

H = — ZiAiajdi + huyi&c + ha>2 (cdi + c^a±j . (15) 

Here, A x = lu d — toi ~ cum, and a 2 = T^/N 2 ~, where N 2 is the average number of photons 
in cavity mode 2. Alternatively, if mode 1 is driven at a frequency wd ~ wi, we get the 
beam-splitter Hamiltonian 

H = — hA 2 d 2 d 2 + huyic^c + hot\ (ca\ + c^a 2 ) , (16) 

in the frame rotating at wd with A 2 = wd — ^2 ~ -wm and a\ = T^/W[ with Ni the number 
of photons in cavity mode 1. It is well known that these Hamiltonians can be realized in the 
resolved sideband regime in ordinary single-mode cavity optomechanics by driving at either 
the blue or red sideband frequency. The difference in this case is that the drive frequencies 
are on resonance with other cavity modes, similar to the situation in Ref. Il2~0ll . 

Finally, we mention that one can also realize a beam-splitter Hamiltonian between the 
two cavity modes and d 2 by driving the mechanical oscillator. The mechanical motion can 
be induced either mechanically or by strongly driving two additional optical modes that are 
also separated by the mechanical frequency. This could allow for transfer of an optical signal 
from one frequency to another. 



5. Estimation of the coupling constant 

In this section, we provide an estimate for the coupling strength between the dilational motion 
and the cavity modes. The size of the zero point fluctuations is independent of the membrane 
thickness for b ^> d and is given by i] zp{ « 2^Jh/ (npb 2 v s ). The speed of sound in Si 3 N 4 is 
v s = 9900 m/s and the density is p = 3.44 g/cm 3 . Assuming that the membrane side length 
is b = 0.1 mm, we get r] zpi = 2.0 x 10~ 17 m. A smaller b will reduce the membrane mass, 
which gives a larger zero point motion and hence a larger optomechanical coupling. This is 
yet another argument for using a small membrane. 

We will imagine that we are using laser light with a wavelength of A = 1064 nm. 
For membrane thicknesses of do = 100 nm and d = 1 pm, the mechanical frequency is 
ojm = 27r x 50 GHz and ojm = 2tt x 5 GHz, respectively. This is relatively large compared 
to many other optomechanical realizations, although some experiments feature mechanical 
frequencies in the GHz regime ll28l |29l [30l . The large resonance frequency means that at 
dilution refrigeration temperatures, the thermal phonon occupation number 



exp(hu M /kBT) - 1 

will be vanishingly small (n th < 10~ 5 even for com = 2tt x 5 GHz), such that the mechanical 
oscillator will practically be in the quantum ground state without the need for additional laser 
cooling. 
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In Section 4.1[ we saw that the optomechanical coupling strength is determined by the 
rate g/2-n = crj zp f/(XL). For a cavity of length L = 1 cm, this becomes 0.7 Hz. This is 
comparable to the coupling strength to the flexural modes of the membrane in experiments 
performed with longer cavities flU. Naively, this suggests that optomechanical effects are 
observable also with the dilational mode. The mechanical linewidth 7 is however expected 
to be much larger for the high frequency dilational mode than for the low frequency flexural 
modes. What we need to check is how the effective coupling rate a ~ g\/N compares to the 
cavity decay rate k and the mechanical oscillator decay rate 7. Here, N is the average number 
of photons in the driven cavity mode. A figure of merit is the cooperativity C = Aa 2 /(ivy). 
Introducing the mechanical quality factor Q = uJu/l an d the finesse T = ctt/(Lk), we arrive 
at 

C-l^Y^NTQ (18) 
V A / v s L 

as an order of magnitude estimate for the cooperativity. 

Inserting the expression for r] zp f, and expressing the average photon number iV = 
4P in A/ (he k) in terms of the laser power P in , our estimate for the cooperativity becomes 

128 d P iQ J^Q 

IT 6 A pb A CVj 

The only unknown in this estimate is the mechanical Q. This can be experimentally 
determined by for example driving the mechanical motion with two laser beams separated by 
the mechanical frequency and analyzing the response of the oscillator. Assuming d — 1 /mi, 
P in = 1.5 mW, and T = 10 5 , we get C/Q ~ 5.8 x 1CT 5 . This means that for a mechanical 
Q > 10 4 , the cooperativity can be on the order of 1, such that optomechanical effects should 
be observable. 



To realize the model (12), we want the frequency difference between the two optical 
modes to equal the mechanical frequency, i.e. d> i+1 = %. In the absence of the membrane, 
the free spectral range of the cavity is Si = ttc/L. Thus, the length of the optical cavity will 
be determined by the requirement L/d ~ c/v s . Let us assume a membrane thickness of 
do = 1 /im, which means that the cavity will need to be approximately L = 3 cm. In practice, 
since the difference 5i depends on the membrane position xm (see Fig. [2]), one can move 
the membrane to the appropriate position where 5i exactly matches %. Note that with this 
frequency matching, the resolved sideband parameter um/k is roughly equal to the finesse T 
of the optical cavity. 

With the assumptions above, the cavity linewidth is k = 2ix x 50 kHz and a ~ 
2n x 60 kHz. For a mechanical quality factor of Q ~ 10 4 , we have 7 = 2rr x 0.5 MHz. 
Thus, even for modest cooperativity, since the effective coupling a exceeds k, one could be 
able to observe strong coupling phenomena such as normal mode splitting of the mechanical 
resonance [[3T1I321 . Significant resolved sideband cooling of the mechanical motion requires 
C 1, which, with the numbers we have assumed here, is only possible if Q 10 4 . 
However, as noted earlier, the thermal occupation will be vanishingly small at millikelvin 
temperatures. 
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It is difficult to give an estimate for what the mechanical quality factor Q of the dilational 
mode will be. It is however worth mentioning that extremely large Q factors exceeding 10 8 
have recently been realized in similar structures, by making one of the surfaces convex, so as 
to provide confinement of the acoustic wave ll33~l . One might imagine that a similar technique 
can be used to engineer large quality factors in the thin membranes we consider here, although 
one should note that the confinement reported in Ref. [f33l is best for higher overtones of the 
fundamental mechanical resonance. For such overtones, where the mechanical frequency is 
com ~ qnvs/do with q > 1 being an integer, the photoelastic term (H 2 ) will give a negligible 
contribution to the optomechanical coupling. The surface term (Hi) could however give rise 
to an optomechanical coupling for odd q. The coupling will still be reduced by a factor 1/ ' ^Jq 
compared to the fundamental mechanical mode, due to smaller zero point fluctuations. To 
confine the fundamental q = 1 mode and avoid propagation losses, it would be preferable to 
use a membrane whose thickness is sharply reduced at the edges. 

Another issue one might worry about is optical loss in the membrane. Note however 
that membrane loss could not be observed in experiments with 50 nm thick membranes and 
a cavity finesse of 5 x 10 4 lfl4ll . While this does not exclude the possibility that it could be 
a problem at a thickness of 1 fim, it is worth pointing out that some of the absorption and 
scattering in the membrane is almost certainly due to surface roughness and contamination, 
which should not scale at all with thickness. 



6. Quantum limited optical amplification 

We will now show how the optomechanical coupling to the dilational mode of the membrane 
can be used for quantum limited optical amplification. We consider the situation where a 
strong laser at frequency wd drives the upper of two cavity modes separated by the mechanical 



frequency, which gives rise to the effective Hamiltonian ( 15 ). The detuning between the drive 
and the lower cavity resonance will be denoted Ai = w D — 0Ji ~ wm. In addition to the 
strong drive, we imagine that a narrow bandwidth signal centered on lo\ is sent to the cavity, 
represented by the operator e~ 1Wl 'ai iin (t). The bandwidth B should be smaller than the decay 
rates k, 7, which for the parameters in Section [5] means roughly B < 50 kHz. The decay rate 
of the mirror through which the laser couples to the cavity modes is denoted by n ex . We will 
account for other decay channels, e.g. leakage through the other mirror or absorptive losses in 
the mirrors, whose respective decay rates are Kj with j integer, such that k = k cx + Y.j K j- In 
a frame rotating at uj\, input-output theory j34l|35l gives the quantum Langevin equations 



K 



at = - -Si - ia 2 e 1 ll c ] + y^ai.in + V^i i 



»3 



\ + i^ M ) c - i« 2 e" iAl * a\ + ^7 C • (20) 

Here, the operator ((t) represents the noise from the mechanical bath. For a large mechanical 
Q, we can use a white noise model with (C(*)C(*0) = n th<5(t — t') and (CWC^O) = 



(n t h + 1) 5(t — t'), where n th is the thermal phonon number defined in Eq. ( 17 ). The operators 
£j represent the vacuum noise of the electromagnetic field and obey = an ^ 
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(ij(t)C k (if)) = Sjk6(t — t'). The output field of cavity mode 1 is ai >out (t) = ^/K^ai(t) — 
a-i,in(t) in a frame rotating at 

We define the Fourier transform as = Jf^ dt e lwt f^\t). For frequencies 

u<k,7, the output field is 

A r , (2r - 1 + C)a X iM + 2^r £j J^Jk^] - iVC&[u>] 

ai,outH = - i-c — — ' 

where r = k cx /k and C = 4a|/(/«7) is the cooperativity. The result is not valid for C > 1, 
where nonlinear terms must be taken into account and self-induced mechanical oscillations 
occur [|36l . In the regime 1 — r < C < 1, the interaction with the cavity leads to a 
phase-preserving amplified signal. Such an optomechanical amplifier has been realized in 
the microwave regime ll38l . although the quantum limit of added noise was not reached and 
the bandwidth was only on the order of kHz. The high-frequency dilational mode discussed 
here has the advantages that the thermal occupation n t h can be made vanishingly small and 
the bandwidth can be significantly larger, due to a presumably larger mechanical linewidth. 
Defining the arbitrary quadrature X e ^ n = (e~ l6 ai )in + e l9 a[ )in )/\/2, we get 

(X e 2 out )=G((X, 2 in ) + (AX) 2 ) , (22) 

where the amplifier gain is 



V 1-C 




and the added noise is 

< A - Y > 2 -kk-hhG- 1 ) i l+ is J (24) 

when nth <C 1. If other types of cavity decay can be neglected, i.e. if k cx ps k, we have 
(AX) 2 = (1 — l/G)/2 and the amplifier is quantum limited ll39l . In the general case and for 
large gain G ^> 1, the added noise becomes (AX) 2 = 1/r — 1/2 > 1/2. This means that 
in practice, deviations from the quantum limit will be due to amplification of vacuum noise 
entering the cavity through other decay channels, and not due to the mechanical oscillator. 

Note also that for cooperativities C < 1 — r, the interaction with the cavity actually leads 
to absorption, in the sense that the gain G is less than unity. In this regime, the magnitude of 
the reflected signal can even become smaller than for no optomechanical coupling (C = 0) 
ll28ll . and a portion of the signal is lost to the other optical decay channels. This has been 
observed in an optomechanical system [|28l , and is analogous to electromagnetically induced 
absorption in atomic gases ll37l . 



7. Correlations between cavity modes 

In this section, we show how the dilational mode of the membrane can be used to create two- 
mode squeezing between different cavity modes. We also show how this system can lead to an 
explicit observation of quantum back-action of the optical field on the mechanical oscillator 
POll . also known as radiation pressure shot noise. 
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Let us now consider four optical modes hi, i — 1, 2, 3, 4, where a 2 (a 4 ) and a x (a 3 ) differ 
in frequency by the mechanical oscillator frequency w M such that scattering between these 
modes can occur through phonon creation or annihilation. In other words, we consider the 
Hamiltonian 

4 

H — hujia\di + huuc^c + hT (c a 2 ai + c*a\a 2 + c 0.40,3 + c^a^a^j . 

Note that we do not require a 2 and a 3 to be adjacent modes, although they may be. It is 
also not necessary that oj\,U)2 < ^3,^4- We imagine that both modes 2 and 3 are strongly 
driven at frequencies u Dj2 ~ cu 2 and w Di3 ~ co 3 , respectively. This means that the mechanical 



oscillator's interaction with mode 1 is of the parametric amplifier type (Eq. (15)), whereas 
the interaction with mode 4 is of the beam splitter type (Eq. ([T6])). This setup can lead to 
correlations between the modes 1 and 4. The reason is that if the mechanical dissipation 
is sufficiently small, one can adjust the drive strengths such that almost every phonon that 
is created by a downconversion of a photon from mode 2 to mode 1 will be destroyed by 
upconverting a photon from mode 3 to mode 4. In other words, the creation of a photon in 
mode 1 is almost always followed by a creation of a photon in mode 4. We now study how 
this can lead to measurable quantum correlations. 

7.1. Equations of motion and self-consistency 

We will assume that (X>d, 2 — — % C and uj a — aj D3 — oj m <S K , 7> which should not be 
a problem for k, 7 in the kHz-MHz range. We transform the cavity and oscillator operators to 
frames rotating at their respective resonance frequencies, i.e. Oj — > e~ lulit ai and c — > e~ lbJMt c. 
Moreover, we can combine the vacuum noise entering cavity mode i from all decay channels j 
in the operator ^(t) = K~ 1 / 2 ( v /fi^^ i ex + J2j This leads to the equations of motion 



cii = - -01 - ia 2 c t + (25) 
a 4 = — -04 — ia 3 c + y/K £ 4 



- |c - i [a 2 a\ + a 3 a 4 ) + ^7 C 



for cavity modes 1 and 4 and the mechanical oscillator. The average amplitudes a 2 = Ta 2 
and a 3 = Ta 3 with = (oj) have been assumed real and positive without loss of generality. 
These amplitudes should be determined self-consistently, which leads to the equations 



_ _ Ta 2 

2,0 (« + 7) (/«7/4 + a| - al) 
_ _ _ Ta 3 

«3 — G-3,0 — 



. . (2/€ + 7)a 3 

7(n th + l) + 



k{k + 7)/2 + a\ — a 2 



(26) 



(« + 7) («7/4 + al - a|) 



(2«: + 7)a 2 



k(k + 7)/2 + a 3 — a\ 



where a^o = 2(n cx Pi /huiK 2 ) 1 ^ 2 is the amplitude in the absence of optomechanical coupling, 
and P 2 (P 3 ) is the power of the beam driving cavity mode 2 (3). To avoid the mechanical 
oscillator becoming unstable 11361 , we focus on the case where a 3 > o 2 . Also, we limit the 
discussion to the weak coupling limit r<7,K and on not much greater than k, relevant to the 



Optomechanics with a high-frequency dilational mode in thin dielectric membranes 



16 



experiment we propose. In that case, the corrections to a^o on the right-hand side of Eqs. (26) 
are small. Thus, it is a good approximation to replace aij in these terms by a ii0 = Tdi . 

While the correction to a 2 and a 3 induced by the optomechanical coupling is negligibly 



small in the regime we study, we will see that Eqs. (26 ) provide insight to the discussion about 
two-mode squeezing in the next section. 

7.2. Two-mode squeezing 

The fact that a photon creation in mode 4 almost always follows a photon creation in mode 1 
can give rise to peculiar photon statistics, similar to what was studied in Ref. ||4~T|. Here we 
study how it leads to two-mode squeezing. We define the internal cavity quadratures 

X dl (t) = ^(a 1 (t)e^+a\(t)e^) , (27) 

From these, we can define a general intermode quadrature 

X 6l>eA {t) = ~(x ei {t) + Xe A {t)) . (28) 

We note that in the absence of the mechanical oscillator, the vacuum fluctuations lead to 
the standard quantum limit (Xg) = (Xg) = (Xg e ) = 1/4 for the variance in these 
quadratures. 

Taking into account the coupling to the mechanical oscillator leads to squeezing in the 
quadrature Xg lt g 4 for appropriately chosen phases Q\ and 4 . Relative to the choice of defining 
a 2 and a 3 real, the optimal choice of phases is d\ + # 4 = modulo 2ir. This results in the 
variance 

= il 1 + (^ + 7)( ( r7/4 + i|-a|) (29) 



7 (("3 - a> 2 )n th - a 2 J 



(2k + 7)a2«3(«3 - Oi 2 ) 
k(k + j)/2 + a\ — a\ 



By using Eqs. (26), we can rewrite this result as 



(Xg ± _ 01 ) = j{l + (o 3 - 02) (03,0 - o 3 ) - (a 2)0 - a 2 ) }, (30) 

which shows that the two-mode squeezing occurs when the average amplitude deviation 
«i,o — a>i is larger for mode 2 than mode 3. From the discussion in Section |7TTj it is clear 



that we can replace a, in Eq. (29 ) by a i)0 in the regime we are interested in where r <C kj, 
with a-sfl > a 2 fl and a it0 not much greater than k. 

The variance (29) is shown in Figure [4] as a function of the ratio between the drives, 



a 2,o/ Q; 3,o- The vertical axis is normalized to the shot noise level of the unsqueezed state. We 
use the parameters estimated from Section[5j choosing n/a^o = 5/6 and nth = 0. We plot for 
three values of the mechanical linewidth, 7 = 10/t, k, and k/10, corresponding to Q = 10 4 , 
10 5 , and 10 6 . We observe that the variance goes below the shot noise level, indicating two- 
mode squeezing. On the other hand, the individual variances (X| ), (Xg ) > 1/4, i.e. there is 
no intramode squeezing. 
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0.4 0.5 0.6 

a2,o/a3,0 



Figure 4. The variance of the quadrature Xg 1 normalized to shot noise level, showing two- 
mode squeezing between modes 1 and 4. On the horizontal axis is drive strength at the blue 
sideband of mode 1 (0:2,0) relative to drive strength at the red sideband of mode 4 (0:3,0)- We 
have assumed K/03,0 = 5/6, n t h = 0, and plot for three levels of the mechanical linewidth, 
7 = k/10, k, and 10k, corresponding to Q = 10 4 , 10 5 , and 10 6 . 



The experimentally accessible quantities are the output quadratures 

X ei ,out(t) = \ (ai,out(t)e- iei + a{ out (t)e^) (31) 
^ 4l out0) = \ (a4,out(t)e- i94 + a{ out {t)e m ) , (32) 

where the output fields are a iiOUt (t) = y/K^a^t) — &,ex(£) an d ii,ex(t) is the incoming vacuum 
noise at the input/output port where the quadratures are measured. These quadratures can 
be measured via homodyne photodetection, yielding access to the (two-mode) squeezing 
spectrum B21B5I 

roo 

S sq [u] = 1 + 8 dr (: X eii e 4iOUt (0)X eiie4iOUt (r) :) cos(wr) (33) 
Jo 

where X Sl t g 4 )OUt = (Xg uont + Xg 4 out ) / and the colons indicate normal and time ordering. 
The spectrum is normalized in such a way that unity refers to shot noise level. We have 
neglected the fact that the two local oscillators used to measure Xg ltOUt and Xg iiOUt may have 
different strengths, and assumed a detector efficiency of unity. However, these are purely 
technical issues that can straightforwardly be taken into account [|43l . Exploiting the fact 
that £i )CX (t) is vacuum noise and that it cannot be correlated with the intracavity operators 
ai(t'), a,i(t') at earlier times t' < t due to causality, the squeezing spectrum can be expressed 
in terms of the internal operators: 

S sq [u] = 1 + 8 n ex / dr (: X ei . e4 (O)X 01 ,, 4 (r) :) cos(u;r) . (34) 
Jo 

Focusing again on the optimally squeezed quadrature (#1 + # 4 = modulo 2n), we find 

S sq [u) = 1 + k cx (a 3>0 - « 2 ,o) IxcMHxmMI 2 (35) 
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X {t[( q; 3,0 - a> 2 fi)n th - «2,o] - ^02,0^3,0 ("3,0 ~ «2,o) |Xc[^]| 2 } . 



We have defined the susceptibilites 

1 



and 



XmM 



k/2 — ioj 



7/2 + (o§ - al)xcM ' 



(36) 



(37) 



The validity of the expression (35) is also limited to the weak-coupling regime discussed 
above. 
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Figure 5. The two-mode squeezing spectrum S^a;]. On the vertical axis is £12,0/0:3, o> 
i.e. drive strength at the blue sideband of mode 1 relative to drive strength at the red sideband 
of mode 4. On the horizontal axis is frequency ui in units of k. We have chosen k/ 013.0 = 5/6 
and n th = 0. Top left: j/k = 0.1 (Q = 10 6 ). Top right: j/k = 1 (Q = 10 5 ). Bottom left: 
7/k = 10 (Q = 10 4 ). Bottom right: j/k = 100 (Q = 10 3 ). 



The output squeezing spectrum is plotted in Figure [5] as a function of frequency to and 
the ratio between the drives, «2,o/ a 3,o- We have used the same parameters as in Fig. [4| and 
plot the spectrum for Q = 10 3 , 10 4 , 10 5 , and 10 6 . We see that the spectrum goes well below 
shot noise level and that nearly maximal squeezing can be achieved at low frequencies to < k 
when the mechanical linewidth 7 < 0:3,0, k and the ratio of the two drives is appropriately 
chosen. Even for a modest mechanical Q = 10 3 , the squeezing spectrum goes several percent 
below the shot noise level. 

It is clear from Eqs. ( [29] ) and (35) that the two-mode squeezing will diminish as the 
temperature and hence the equilibrium phonon number n th increases. For extremely large 
mechanical Q, the squeezing could nevertheless be observable at higher temperatures. For 
example, to make the variance (29) go 1% below the shot noise level at room temperature 
would require Q > 3 x 10 7 with the parameters we have assumed. 
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7.3. Radiation pressure shot noise 

Observation of the squeezing discussed above implies that the mechanical oscillator is 
susceptible to radiation pressure shot noise. We now discuss another way of using the two 
quadratures Xg lfiUt and Xe 4:OUt to explicitly detect these correlations between the mechanical 
oscillator motion and photon shot noise. The method we propose is similar to the detection 
schemes studied in Refs. P4l |45ll . However, the experimental setup we have studied here 
seems particularly appropriate for this kind of measurement, warranting a thorough analysis. 



As in the previous section, we focus on the model in Eq. (25 ) and imagine that modes 2 
and 3 are strongly driven. Homodyne detection of Xg ljOUi (t) and X$ i}0n t(t) gives access to the 
cross-correlation spectrum 

SM = \T dre iwT ({l eij0ut (r) , X e4 , out (0)}> , (38) 

Z J — CO 

where {•, •} denotes the anticommutator. This spectrum can be expressed in terms 
of two types of correlation functions. Some correlation functions are of the type 
($(t)c(0)), i.e. mechanical oscillator autocorrelation functions. The others are of the type 
(£j(r)c(0)), which are quantum optomechanical correlations between the vacuum noise of 
the electromagnetic field and the mechanical oscillator. The latter type is nonzero when the 
mechanical oscillator is susceptible to photon shot noise. In other words, these correlations 
functions are nonzero due to radiation pressure shot noise. 
The cross-correlation spectrum becomes 

S cr [u] = -K cx a 2 ,oa3,o|xcM| 2 |xM[w]| 2 cos(6 l i + d 2 ){R CT [uj} + LT C rM) , (39) 



where 



and 



R cr [u] = 7 (2n th + 1 ) + k (a 2 . + a 2 ) | X c M | 2 (40) 



I CI [u] = -2a; [l + (a 2 >0 - a 2 2fi ) |xcM| 2 ] . (41) 

Note that the real part of S ct [uj] is even in frequency, whereas the imaginary part is odd. These 
functions are displayed in Figure [6] for the case of equal driving strengths, a 2i0 = a 3i0 . The 
real part comes from the oscillator autocorrelation functions, whereas the imaginary part is 
solely due to the correlations between the mechanical oscillator and the vacuum noise. Thus, 
measurement of an odd frequency imaginary part of the correlation function S ci [uj] is clear 
evidence of radiation pressure shot noise. 

In the calculation of S cv [u] , we have neglected technical laser noise, which can give rise 
to a similar asymmetry as the vacuum noise if it is not negligible at frequencies around wm 
B3I . However, for a mechanical frequency in the GHz range, neglect of the laser noise is a 
very realistic assumption. 

Finally, it is worth noting that the antisymmetric behaviour of the imaginary part I ct [uj] 
can in principle be measured at high temperatures. However, as the magnitude of the real 
part is proportional to temperature, control of phases so as to correctly identify the real and 
imaginary parts would have to be very accurate. Challenges of this type were discussed in 
more detail in Ref. P51l . 
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Figure 6. The real (left) and imaginary (right) parts of the cross-correlation spectrum 5 cr M- 
We have used n/a^.o = 5/6, a 2 , 0/^3,0 = 1. j/k = 1, and n t h = 0. The asymmetric 
imaginary part is a signature of radiation pressure shot noise. 



8. Concluding remarks 

We have studied the optomechanical coupling between the thickness fluctuations of a thin 
dielectric membrane and an optical cavity field. For an appropriate choice of membrane 
thickness versus cavity length, the model reduces to one where the creation (annihilation) 
of a phonon causes scattering from an optical mode with higher (lower) frequency to one with 
lower (higher) frequency. We have derived the coupling constant and showed that realistic 
parameters suggest that optomechanical effects can be observable. We argued that this system 
can be used to amplify a narrow bandwidth optical signal with quantum limited added noise. 
Finally, we showed how one can create quantum correlations between two optical cavity 
modes at different frequencies. We calculated a two-mode squeezing spectrum and found 
considerable squeezing for realistic parameters. A method for clearly detecting radiation 
pressure shot noise was also presented. 

The physical system we discussed is already being actively studied experimentally 
in several labs, such that very small changes are necessary to attempt to observe these 
phenomena. A primary challenge is to detect the dilational motion, either through its thermal 
noise or its response to driving, and to determine the mechanical quality factor. If the latter 
is small, it would be an important technical challenge to find ways to increase it. With a 
sufficiently high quality factor, the possibility of observing quantum optomechanical effects 
with this high-frequency dilational mode appears feasible. 
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